Deterministic Weak Localization in Periodic Structures 
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The weak localization is found for perfect periodic structures exhibiting deterministic classical 
diffusion. In particular, the velocity autocorrelation function develops a universal quantum power 
law decay at 4 times Ehrenfest time, following the classical stretched-exponential type decay. Such 
deterministic weak localization is robust against weak enough randomness (e.g., quantum impuri- 
ties). In the ID and 2D cases, we argue that at the quantum limit states localized in the Bravis cell 
are turned into Bloch states by quantum tunnelling. 
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The remarkable phenomenon of strong localization 
(SL) ll| nowadays has been fully appreciated in aperi- 
odic and periodic 0] disordered electronic systems and 
disordered dielectric materials ]3| . The weak localization 
(WL), as the precursor of SL [4j has been confirmed by 
tremendous experiments in the last two decades, e.g., the 
negative magnetoresistance , the coherent backscatter- 
ing of light p], etc.. All these systems are strongly scat- 
tered, i.e., A > a, where A is the wavelength involved 
and a is the typical scale over which the potential varies 
substantially. Yet, in nineties unprecedented degree of 
control reached in experiments with mesoscopic quan- 
tum dots allows to investigate the semiclassical region: 
A <C a , where quantum interference is suppressed, and 
classical chaotic motion prevails. Among central issues is 
the classical-to-quantum crossover. A long time ago Q 
it was established that appears the so-called Ehrenfest 
time marking the proliferation of quantum interference. 
However, the quantitative crossover had not been found 
until recently |9(, in the context of quantum corrections 
to the Drude conductivity of disordered ballistic systems. 

It is well known P, Q that in disordered time-reversal 
symmetric systems, WL originates from the quantum 
interference between two counter-propagating trajecto- 
ries, and the diffusive nature of trajectories arises from 
the randomness, regardless of classical or quantum po- 
tentials. In contrast, the deterministic diffusive motion 
may take place on a classical periodic potential, where 
the extended Bloch state develops at the quantum limit. 
Transport properties in such system remains unexplored. 
On the other hand, the so-called periodic Lorentz gas, 
as a pro toty pe has been well understood mathematically 
pLOL rLlL Il2fl , where the interaction with the potential is 
simplified as "mirror reflection" . Essent ial prop erties are 
encapsulated in a series of theorems [IollllLll2| . To sum- 
marize, (i) the flow is K mixing |ic|: (ii) for piecewise 
Holder continuous (PHC) functions the correlation decay 
[e.g., the velocity autocorrelation function (VCF)] is fast 
of stretched exponential type, i.e., ~ exp(— £ 7 ) , < 7 < 1 
and (iii) the diffusive dispersion relation is identi- 



fied as a Pollicott-Ruelle resonance with an exact Green- 
Kubo relation for the diffusion constant established ^3] • 
A finite diffusion constant is ensured by a priori, so-called 
finite horizon condition |lfjj . To explore quantum man- 
ifestations of these features naturally becomes a funda- 
mental problem. 

The main finding of the present work is to show that 
such deterministic classical diffusion leads to WL (to dis- 
tinguish from the usual WL in disordered systems, here 
we may term it deterministic WL). In particular, at the 
one-loop level, we find the frequency-dependent quantum 
correction to the diffusion constant to be 

^)^-^rn/A : 1 (!) 

■kTiv J (2tt) z —iuu + DclK.^ 

where v is the density of states, and the renormalization 
factor 



T(uj) — exp AiuttE — 



The Ehrenfest time, say t E is 



t E = A" 
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Here A is the Lyapunov exponent and can be expressed in 
the form of Abramov formula fl^ | . Depending on initial 
conditions, t E may fluctuate, characterized by A2 with 
the characteristic scale StE = A2£,e/A 2 . Results similar 
to Eqs. Q-lEJ were found for disordered ballistic sys- 
tems earlier |9|. However, we emphasize that throughout 
the derivation below no quantum impurities will be in- 
troduced. Therefore, the main problem surrounding the 
result of Ref . |fj , namely the possibility of removing aux- 
iliary quantum impurities is solved. On the other hand, 
for the first time we see that at the quantum limit WL 
leads to extended states rather than SL. 

Alternatively, according to the Green-Kubo formula 
one may view the result above from VCF (v(t) • v(0)) 
[with the average over the primitive cell of the phase 
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space energy shell (for shorthand we will call primi- 
tive cell below.)]. In the semiclassical limit, tg satis- 
fies Tfi <C i_E "C th , where the localization time scale 
~ mdy 1 D c i/h and Tfi is the classical mean free time. 
Ignoring St E , we find that the VCF develops a power 
law decay at t > At E beyond the classical stretched- 
exponential type decay. Particularly, in the 2D case, the 
power law takes the form as — (Air 2 his)" 1 (t — Ats) -1 (see 
Fig.Q). 

<v(t) v(0)> 



expf—f 1 ] 
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FIG. 1: The velocity autocorrelation function of the motion 
in a 2D periodic Lorentz gas (with finite horizon) develops a 
power law decay at t > <Ue ■ 

Remarkably, the standard WL correction acquires a 
dispersion T(uS) , prevailing at u> ~ t E x . In fact, simi- 
lar renormalization factors exist at all the higher order 
loop corrections. Such phenomenon originates from the 
interference nature of the localization. Indeed, two clas- 
sical trajectories pass through (almost) the same posi- 
tion at (almost) the same momentum. Subsequently they 
diverge and eventually take counter-propagating routes. 
If these positions and momenta were strictly identical, 
the probability would be zero. The quantum uncertainty 
makes it possible. It takes, however, a long time ~ t E . 

Having outline the qualitative picture, we turn to de- 
tails of the proof. To be specific, let us focus on the 
2D case. The model consists of an electron of mass 
m with the energy E (namely the Fermi wavelength 
moving in a periodic Lorentz gas. The 
lattice constant is d, while a is referred to the disk ra- 
dius. Technically, to avoid the singularity of the hard 
core potential, the billiard is softened to be the poten- 
tial V(r) around the boundary with the typical width 
b -C W := d — 2a . The Hamiltonian thereby is replaced 
by H -.= p 2 /2m + F(r) . 

The basic tool that we will employ is the gen- 
eralization of the diagrammatical technique de- 
veloped in essentially different context-the kicked 
rotor 13]. We will skip parallel intermediate steps, 



with emphasis on the main difference. Let us start 
from the quantum four-point "density-density" 
correlator, defined as: 2?(r+, r_; r' + , r'_; t, t') :— 
(r + | exp[irlt/h]\r' + ) (r_| cxp[iHt' /h]\r'_}* , where 
r± , r'± . A crucial step is to introduce an artifitial 
"one-step" evolution operator U := exp[i7ii*/7i] , and 
consider the density-density correlator at times of 
multiple i*, i.e., t = nt* ,t' — n't* . Here t* is a time 
scale such that <§; Tfi . Passing to the frequency 
representation-^ , E) , the density-density correlator is 

oo 

2? £ (r + ,r_;rV,r'_;c) = £ <r + |?7"|r' + ) <r_|l>"'|r'_)* 
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Now we consider the classical limit of the exact quan- 
tum density-density correlator. First let us sum over all 
the diagrams such that |r + — r„ | , |r' + — r'_| , |r^ + — rfc_| ~ 
Af <C W , k = 1, 2, • • • [see Fig. [21 (a)] and denote it as 
T>o . To proceed further, we perform Wigner transform 
the quantum density-density correlator with respect to 
r+ — r_ , r' + — r'_ and pass to the (p , p') representation. 
As a result, we find that T>o evolves following the clas- 
sical Perron-Frobenius equation according to (u> will be 
ignored to shorten notations.) 



to, ( iuo — C] 



1 — exp 



(2tt^) 2 6 (r - r') 5 (p - p') 6 IE 



)]}p B0 (r,p;r',p') (5) 
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where C = v ■ V r — V r V(r) ■ V p is the Liouvillian. Note 
that in most region of the phase space the Liouvillian is 
free except near the boundary of the billiard, where the 
potential is involved. According to Eq. JSJ, the modified 
Hamiltonian induces a flow restricted on the energy shell: 
E = p' 2 /2m. We assume that for this flow, the essen- 
tial properties (i), (ii) and (iii) for the periodic Lorentz 
gas hold. To be specific, in the PHC functional space, the 
Fourier transform of the operator {1 — exp[t*(iw — £)]} _1 
admits the diffusive poles: u> — —iD c ik 2 (k the wave 
number taken from the first Brillouin zone of the recip- 
rocal lattice). Moreover, the classical diffusion constant, 
D c i is given by the Green- Kubo formula, i.e., 



D d = 



4 m 2 



E (pCj)-p(o)) 
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with the average over the primitive cell. We point out 
that D c i is finite due to the nature of the finite hori- 
zon, although we leave the problem of calculating Eq. © 
open. Based on this conjecture, despite of the determinis- 
tic nature of Eq. (0) , one may work on distributions with 
the smoothness admitted by the PHC condition. Define 
V E0 (r, p; r', p') := V (r, n; r', n') 5(E - p 2 /(2m)) 8{E - 
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p /2 /(2m)) . Averaging Vq over the primitive cell, and with 
passing to the Fourier representation, eventually we re- 
duce T> into the diffusive propagator (denoted as T> v ) 



V v k;u, = 



(7) 



(a) 



(b) 



(c) 
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D, 



= exp[ it *(o]uu + - l 



: retarded propagator 
: advanced propagator 



FIG. 2: Typical diagrams lead to the Diffuson (a) and the 
Cooperon (b). The difference of two coordinates connected 
by the dot-dashed line denotes is order of Af <§C W . The 
one loop correction is given by (c). The retarded (advanced) 
propagation line is exp[it»a)+]?7 (exp[—it*Lj-]ffi). See Eq. © 
for the definition of *C and C, . 



t"_| , |r' + -r_|, |r fe+ -r (n+1 _ i . ) _| ~ A F < 
n [Fig. [21 (b)] yields another classical 



We note that in fact, to sum over all the diagrams 
such that |r+ ■■' 
W,k = 1,2 

object, so-called Cooperon (denoted as C ). In fact, one 
may Wigner transform the quantum density-density cor- 
relator with respect to r + — r'_ , r' + — r_ and arrive at 
Eq. ©. 

Having identified the classical normal diffusive modes, 
we come to analyze WL. For this purpose, the essen- 
tial step is to show that in the semiclassics, the one-loop 
quantum correction is given by the diagram shown in 
Fig. [21 (c). In the exact quantum operator representa- 
tion, it may be expressed as 



5V = V 



EQ 



(. 



c + c, 
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C» := X^VC E oT n . (9) 



Here V = e iu)t * UU^ . Notice that in Eq. © the operator 
V in (V — 1) is chosen in the way such that it describes 
the free motion. As usual, we then perform the Wigner 
transform for T>eo , *C , C* and V, respectively. Now we 
assume that the initial distribution is PHC. Furthermore, 
the appearance of diffusive poles allows us to make the 
so-called hydrodynamic expansion, u>t* <C torfi <C 1 and 
exp[i*(iw — £,)] - la — £*v • V . Eventually in the PHC 
functional space, the one-loop quantum correction to the 
Diffuson is reduced into 



W?(r,n;r',n') = ^V 
2 m 
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where n(0) := n" . The vertex operator V is an integral 
operator defined as 

dr"dn" f dSridSiii f dSr 2 dSn 2 



V/ (•) := 



2-Kti 



2wh 



x C | )• 

x / (.;r",n";*r 2s <yni) 
where 



Sri „ 5n 2 „ 5^ „ Sn 2 \ 
— — ,n — :r — , -n — 
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(<5ri • 5m + 5r 2 ■ Sn 2 ) 
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(12) 



Above the dot is the shorthand of (r, n; r', n') . The inter- 
ference factor X originates from the phase difference of 
two paths underlying the Cooperon, and technically re- 
sults from the difference of performing Wigner transform 
for the Diffuson and the Cooperon. X can be also ob- 
tained in the Moyal formalism 0]. In the diffusive limit, 
Eq. l|lUfl leads to the standard WL in disordered systems. 
For ballistic systems like periodic Lorentz gases, the two 
legs of the Cooperon must propagate together for 2tE in 
the Lyapunov region in order to develop the universal 
quantum limit of Cooperon develops, which stands for 
the probability of returning to the origin via the random 
walk. Such a crossover picture, encapsulated in the renor- 
malization factor T(u>) was proposed earlier 0. Here the 
essential difference is that, due to the absence of quan- 
tum impurities the motion equation of Cooperon does 
not involve any regularizers. Consequently the Ehrcnfcst 
time is fully determined by the size of the minimal quan- 
(8) turn wave packet, as read out from X , i.e., 5i"i. m j n ~ \ F 



4 



and (friomin ~ Xp/W . Then following the procedure of 
Ref. UM one ends up with Eqs. (01 and JSJ|. 

At the time ~ tp , / acted by V is smeared out over the 
primitive cell surrounding the entries (/ thereby obtained 
is denoted as f v ). Consequently it can be shown that the 
interaction vertex V can be simplified as (denoted as V„) 



(V/) (• 



V = 



Vv/*J (■) :=Vy dr"/„(r,r';r") 
T(lu) f dk 1 



■nhv J (2tt) 2 —iu + D c {k 2 



(13) 



Now let us identify the Diffusons in the right hand side 
of Eq. (|1U[) as deterministic diffusive modes. Technically, 
we may average (r, n) and (r',n') over the primitive cell 
surrounding them. With the application of Eq. ljT3|) . we 
thereby obtain the one-loop quantum correction as 



7 V> 



6V(t,t') = UD cl V J dv" (V^+V* 

[V v (T,T 1 )V„(T 2 ,I f )]\ ri=ra=rlf , (14) 

where we take into account the isotropy of the Bravis 
lattice. Then we perform the Fourier transform and sub- 
stitute Eq. (J7J, as well as the renormalized interaction 
vertex V [see Eq. <|13[l ] into it. As a result, we find that 
the diffusive modes are modified as [— it^uj + t*(D c i + 
SD^k 2 }- 1 with SD(uj) given by Eq. Q. 

An important question is how WL develops into ex- 
tended states in the quantum limit. For qualitative dis- 
cussions, it may be instructive to discuss the quasi-ld 
case. For A_f <C W -C d, each Bravis cell has a well de- 
fined mean level spacing, say A ~ ft 2 /md 2 , with a tun- 
nelling energy, say t between them. Increasing the size 
of the cell n times larger then we rescales the mean level 
spacing as A„ ~ A/n. For t n , we note that the Thou- 
less energy of the rescaled cell is of the order of the total 
tunnelling energy, namely HD c i/(nd) 2 ~ £^/A n , yield- 
ing t n ~ (hDciAn) 1 / 2 /nd . Localization states develop at 
A„ > t n , namely n > n c with n c ~ mD c i/h ~ W/Xf the 
localization length. In contrast to the usual disordered 
system, these localization states (with the same energy) 
are exactly identical because of the periodic feature, and 
thus, serves as Wannier function in the electronic band 
theory. The tunnelling between these states turns them 
into the extended state. In particular, at the quantum 
limit, W ~ Xf , n c ~ 1 , the energy band takes the form 
of Bloch band, i.e., e + t na cos(dk/2ir) . 

In realistic environments randomness (e.g., quantum 
impurities) dooms to exist. Bloch states are unstable 
against randomness. Indeed, in the ID and 2D cases ar- 
bitrarily weak randomness leads to SL at the long time 
(or large size) limit. In contrast, the deterministic WL 
is robust. In fact, any randomness results in a scat- 
tering rate tZ' 1 . As far as it is weak enough, namely 



size of the initial quantum wavepacket due to the loga- 
rithmic accuracy. This conclusion, indeed is consistent 
with an important conjecture made in Ref. Q. That 
is, although auxiliary quantum impurities may be intro- 
duced artificially, the strength of the scattering rate r" 1 
must be adjusted to be ~ {hv)^ 1 in order to mimic quan- 
tum diffractions. Thus, an estimation for the Ehrenfest 
time was established , which is the same as Eq. (j2J . In 
the present work, this essential conjecture is proved. 

To conclude, we find WL for periodic Lorentz gases 
with finite horizon. The result is applicable for any sys- 
tems exhibiting deterministic diffusive motion, regard- 
less whether the underlying potential is ordered or dis- 
ordered. Finally since WL is responsible for coherent 
backscattering |6|, we expect that the result here may 
be confirmed by experiments on coherent backscattering 
of light in state-of-art periodic photonic materials, where 
the semiclassical condition: A <C a may be met. 
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